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Abstract. In situ particle acceleration is probably occur- 
ing in cluster radio haloes. This is suggested by the unifor- 
mity and extent of the haloes, given that spatial diffusion 
is slow and that radiative losses limit particle lifetimes. 
Stochastic acceleration by plasma turbulence is the most 
likely mechanism. Alfven wave turbulence has been sug- 
gested as the means of acceleration, but it is too slow 
to be important in the cluster environment. We propose, 
instead, that acceleration occurs via strong lower-hybrid 
wave turbulence. We find that particle acceleration will be 
effective in clusters if only a small fraction of the cluster 
energy density is in this form. 



1. Introduction 

Diffuse synchrotron emisison from the intercluster plasma 
(ICM) is produced by relativistic electrons moving in the 
cluster magnetic field. The origin of these particles is prob- 
ably easy to explain. They may have been produced ini- 
tially by galaxies undergoing active star formation at early 
epochs (Volk & Atoyan 1999), or they may be injected 
continually into the ICM by current-epoch galactic winds 
and stripping (Petric & Eilek 1999). It has also been sug- 
gested that they can be created by active galaxies in the 
cluster (Feretti et al. 1995), or as secondaries from a rel- 
ativistic baryon population (EnBlin 1999). The particles 
must then diffuse through the ICM; diffusion is a slow pro- 
cess (Eilek 1992, Kirk et al. 1996), but given distributed 
sources the particles should fill a cluster-scale halo (such 
as Coma) in the lifetime of the cluster. 

Whatever their origin, however, their maintenance 
must be explained. The radio-loud electrons lose energy 
by inverse Compton losses on the microwave background, 
with a lifetime ~ 250 Myr. There is also some evidence 
for magnetic fields in excess of 3/xG in parts of some clus- 
ters; such strong fields will shorten the radiative lifetime 
for these electrons. 

Thus, it seems very likely that the electrons are under- 
going in situ reacceleration in the diffuse ICM. How might 
this occur? Shock acceleration is proposed for many astro- 
physical situations. It is not obvious to us, however, that 



shocks are common throughout most clusters. The ICM 
is typically a bit warmer than the galaxies, so only a few 
galaxies will have localized bow shocks. Ongoing merg- 
ers of sub-clusters will, indeed, generate perpheral shocks 
which may lead to localized particle acceleration (Roet- 
tiger 1999, EnBlin et aL1998). However, the diffuse radio 
haloes seen throughout the volume of some clusters prob- 
ably require re-acceleration without shocks. 

It follows that stochasatic, turbulent acceleration must 
be taking place in radio halo clusters. In particular, it 
must be fast enough to offset the radiative losses suffered 
by the radio-loud electrons. Motivated by this question, 
we are studying turbulent particle acceleration in the dif- 
fuse ICM. We are particularly interested in the detailed 
coupling between the turbulent energy density and the 
relativistic particles. It is well-known that acceleration by 
turbulent Alfven waves is slow; in the cluster environment 
it is unlikely to be able to offset radiative losses. We pro- 
pose that some of the turbulent energy resides in nonlinear 
lower hybrid (LH) waves, which are much more effective 
at particle acceleration. In this paper we summarize our 
current work on LH wave acceleration. Details of our ar- 
gument are presented in Eilek, Weatherall & Markovic 
(1999), and also in Weatherall & Eilek (1999). 

2. Turbulent Acceleration, in General 

To begin, we recall the basics of turbulent acceleration 
competing with radiative losses. Turbulent acceleration re- 
sults in a stochasic diffusion in momentum space, while 
radiative losses lead to a direct flow through momen- 
tum space. The governing equation can be written (e.i?., 
Borovsky & Eilek 1986), 



df 1 d f , df 



dt p^ dp 



dp 



Sp'f 



(1) 



The first term on the right describes turbulent acceler- 
ation: Dp is the diffusion coefficient (which can be esti- 
mated as (Ap)^/Ai, where Ap is the mean momentum 
gain in a collision of a particle with a wave or wave packet, 
and At is the mean time between collisions. The second 
term describes radiative losses, if the single-particle loss 
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rate is dp/dt ~ Sp^. Here, S is proportional to the sum 
of the energy density in the magnetic field and the mi- 
crowave background. From this, we can also estimate an 
acceleration time for particles at momentum p: 

r,2 



M^ 



p- 



Dpip) 



(2) 



Analytic solutions of this equation exist for a steady- 
state, closed system with Dp = DoP^ (Borovsky & Eilek). 
They describe a distribution function, 

3- 



f{p) (xp exp 



(3) 



which is peaked at the momentum pc ~ {Do/Sy^^'^^^\ 
where gains balance losses. These solutions are illustrated 
in Figure 1 for r = 3/2 (which may describe Alfven wave 
acceleration), and for r = (which may describe lower 
hybrid wave acceleration). 




Fig. 1. Typical particle distribution function for turbu- 
lent acceleration combined with radiative losses. The dot- 

n3/2 



not damp easily. They have also been suggested as an ac- 
celeration mechanism in other diffuse astrophysical plas- 
mas, such as radio galaxies. They cannot, however, easily 
explain in situ acceleration in the ICM. In this section 
we present a brief overview of the physics of Alfven wave 
acceleration; Eilek & Hughes (1991) contains a fuller dis- 
cussion and references to the original work. 

Alfven waves are long-wavelength, low-frequency 
waves which propagate close to the magnetic field direc- 
tion. Most of the wave energy is in the transverse per- 
turbed magnetic field; only a small fraction is in the elec- 
tric field of the wave (which is capable of doing work on 
the particles). 

Alfven waves interact with particles through the cy- 
clotron resonance, when the (Doppler shifted) wave fre- 
quency is equal to a multiple of the particle's gyrofre- 
quency: u — kuvn + sfl = (where uj is the wave fre- 
quency, A; 1 1 the component of the wave vector along B, fl 
is the particle's gyrofrequency and s an integer. Applied 
to relativistic particles and a typical turbulent spectrum 
of Alfven waves, this is a restrictive condition: a particle 
at momentum p sees only one wavenumber, k ~ mil/p. 
The effect of this is that only a small fraction of the tur- 
bulent Alfven wave energy can be utilized to accelerate a 
given particle. 

The wave-particle interaction has been treated by sev- 
eral authors in the quasi-linear limit. The result from this 
work which wc need here is the diffusion coefficient. This 
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where Uresip) is the turbulent energy at wavelengths res- 
onant with particles at p; note in general that Ures ■^ Ua, 
the total energy density in Alfven waves. 



ted line is the solution for Dp ex p^/ (a likely solution for ^ j^ower Hybrid Turbulence is Fast 



Alfven wave acceleration), while the solid line is for Dp 
independent of p (the probable solution for LHT acceler- 
ation). From Borovsky & Eilek (1986). 

Our discussion up to here has been set in general terms. 
The heart of the turbulent acceleration problem, however, 
lies in determining the diffusion coefficient, Dp{p). Com- 
puting Dp requires understanding the specific physics of 
the interaction between the turbulent plasma waves and 
the relativistic particles. In the rest of this paper we com- 
pare two possible situations, Alfven wave turbulence and 
lower hybrid wave turbulence. 

3. Alfvenic Turbulence is Slow 

Alfven waves are commonly invoked for wave-particle 
trapping and acceleration in the ISM; they are easily gen- 
erated by macroscopic (hydrodynamic) processes, and do 



Lower hybrid turbulence has not often been proposed in 
astrophysical settings, but is well known to exist in the 
terrestrial ionosphere where it is associated with particle 
energization (e.^., Kintner 1992). In this section we dis- 
cuss their effect on relativistic particle acceleration, sum- 
marizing our more detailed argument in Eilek, Weatherall 
& Markovic (1999). We briefiy consider the source of the 
waves in the next section (also Weatherall & Eilek 1999). 
Lower hybrid (LH) waves are short-wavelength, mid- 
frequency waves which propagate across B. They have 
a characteristic frequency, ujlh^ equal to the geometric 
mean of the thermal ion and electron gyrofrequencies. 
Most of the wave energy is in the longitudinal electric field 
of the wave (and thus is available for particle acceleration). 
When LH turbulence becomes strong, it develops intense, 
localized wave packets: field-aligned filamentary regions 
of high electric field. Waves inside these wave packets are 
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undergoing "collapse" - rapid growth of the wave ampli- 
tude accompanied with spatial contraction of the packet. 
The collapse preserves the potential drop across the packet 
{e.g., Shapiro et al. 1993, Melatos & Robinson 1996). 

Particle acceleration occurs through transit time 
damping, which happens when a particle crosses the 
packet in a time equal to the wave period. We expect 
that relativistic particles at all energies will interact with 
all wave packets, once the packet size satisfies 1\\u}lh ^ c. 
Thus, all of the wave packet energy is available for accel- 
eration; and relativistic particles will be selectively accel- 
erated before the nonrelativistic population. We note that 
all particles with w ~ c will interact with the wave packets 
once they collapse to this size; this is quite different from 
the Alfven wave resonance condition discussed above. 




Fig. 2. Localized wave packets in strong Langmuir turbu- 
lence, shown in contours of the electric field amplitued; 
from Weatherall et al. (1983). Strong LH turbulence re- 
sults in similarly localized, but more elongated, regions of 
strong electric field. 



We are in the process of calculating detailed particle 
acceleration efficiencies and diffusion {Dp) coefficients in 
strong LH turbulence. As a preliminary step, we have used 
basic scaling arguments to estimate what will happen. A 
particle gains, on average, Ap ~ 2Trr]eEo/ujLH when tran- 
siting the wave packet, if Eo is the packet electric field. We 
are currently carrying out numerical simulations to deter- 
mine the efficiency ij, as a function of particle energy. For 
our present estimates, we take the rj ^ 1% based on previ- 
ous work for non-relativistic particles (Melatos & Robin- 
son 1996). We scale our estimate to the energy density of 
the turbulence, Ulh, use Eo to find the number density 
of such packets, and from this derive the typical time it 
takes a relativistic particle to move between packets. 



From these arguments, we estimate the diffusion coef- 
ficient to be 



Dp = Dlh 



8W 2 



e ctlhUlh 



(5) 



where r <^ 1 is the aspect ratio of the wave packet (trans- 
verse length over parallel length). We note that Dp is 
nearly independent of particle energy (although the factor 
77 may be a function of particle energy). It is instructive 
to compare this to the Alfven wave diffusion coefficient. 
Writing the particle energy as 7 = p/rrieC, we find 
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This demonstrates that, even in the unlikely case where 
Ures ^ Ulh^ inspection of (6) shows that Dlh ^ Da- 
Lower hybrid wave acceleration is thus a much more ef- 
fective means of accelerating relativistic particles than is 
acceleration by Alfven waves. 

5. Acceleration in the Cluster Environment 

To discuss particle acceleration in the diffuse cluster 
plasma, we pick numbers for a "typical cluster" : B ~ l^G; 
n ~ 10~'^cm~'^; T ~ lO^K (other parameters are given in 
Weatherall & Eilek 1999). This makes the cluster plasma 
weakly magnetized, with the plasma frequency above the 
cyclotron frequency, and the sound speed above the Alfven 
speed). We use these to describe acceleration by Alfven 
and by LH waves. 

For Alfven wave acceleration, we must choose a wave 
spectrum. We assume a Kraichnan spectrum, W{k) oc 
fc~'^", with a forward cascade from the driving scale 
\o — 10 kpc. This wave spectrum gives Uresip) ex p^", 
and the diffusion coefficient can be written Da — DoP^^^. 
We scale the total wave energy density, Ua — J W{k)dk 
to the total magnetic field energy density, Ub = iJ^/Svr.n 
When the steady state has been reached, Alfven turbu- 
lence supports a peak particle energy 



7.^ 9800 (^ 
and takes 



2/3 



Tacc ^ 2.7 X 10« ' ^"^ 



Ua 



2/3 



yrs 



(7) 



(8) 



to reach this energy. Thus, if Ua — Ub, we see that 
Alfvenic turbulence can just manage to account for the 

^ Recall that Ures is the fraction of the Alfven wave energy in 
waves resonant with 7, where Ulh is the total energy density 
in lower hybrid waves 

^ The overall magnetic field strength is probably the best 
observational estimator of the turbulent energy. This scaling 
allows us to consider the fraction of that energy going into 
Alfven or LH waves. 
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synchrotron emission we see. (For comparison, a particle 
at 7 = 10'' radiates at 425 MHz in a /zG field.) 

For lower hybrid wave acceleration, we do not need to 
specify the turbulent wave spectrum, but only the energy 
density in waves. We find the steady spectrum peaks at 



7c 



"-»'(^) 



1/3 



(9) 



which is well above the radio range if Ulh — Ub- This 
peak is reached in 



,(7,) ^ 1.7 X 10 



7/ Ub 



2/3 



yrs 



(10) 



\Ulh, 

while the acceleration time for radio-loud electrons (with 
7r ^ 10*, say) can be considerably shorter: 



= (7r) 



2.1 X 10' 



6 Ub 
Ulh 



sec 



(11) 



Thus, even a very low level of LH turbulent energy, Ulh "^ 
Ub , can still easily account for the radio- loud electrons in 
the cluster. 

Our numerical comparison here recovers the result of 
equation (6): lower hybrid waves are much more effective 
at particle acceleration, for a given level of wave energy, 
than arc Alfven waves. 

6. Connections to Radio Haloes? 

To summarize: we propose that particle acceleration in the 
ICM proceeds via strong lower hybrid wave turbulence. 
We have shown that LHT acceleration is much faster than 
Alfven wave acceleration, given similar levels of turbulent 
energy density in the two wave forms. 

How are LH waves generated? They are easily gen- 
erated by streaming instabilities {e.g., Galeev, Malkov & 
Volk, 1995). As is often the case, the particle beams nec- 
essary for streaming instabilities are not obviously gen- 
erated in diffuse astrophysical plasmas such as the ICM, 
and one must look to hydrodynamic generation. This is 
also possible, as LH waves can be generated by nonlinear 
mode coupling between lower frequency MHD waves. Be- 
cause Alfven waves are easily generated and not strongly 
damped in the cluster setting, we are investigating them 
as a source of LH turbulence. In Weatherall & Eilek (1999) 
we present preliminary numerical work showing that large 
amplitude Alfven waves do, indeed, generate lower hybrid 
waves (through parametric excitation), and that the LH 
waves do undergo collapse to produce localized, intense 
wave packets. We thus suspect that normal hydrodynamic 
processes in the cluster will generate some level of LH 
waves. (Figure 3 shows a cartoon of energy flow in a pos- 
sible turbulent cascade). Our argument above shows that 
only a modest level of LH turbulence can be very effec- 
tive at accelerating relativistic particles, so that we can 
afford inefficiency in the coupling between the large-scale 
turbulent drivers and the small-scale LH wave turbulence. 
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Fig. 3. An illustration of the possible energy flow from 
large-scale hydrodynamic drivers, through Alfven waves 
undergoing a forward cascade, to wave dissipation. Alfven 
waves will undergo modest dissipation when they reach 
wavenumbers resonant with the relativistic particles, and 
they will also undergo dissipation through nonlinear gen- 
eration of LH waves. In this cartoon we have suggested 
that LH generation ends the cascade; however that re- 
mains to be verified with numerical simulations presently 
underway. 



What about the distribution function (DF) of particle 
energies? We noted above that any turbulent acceleration 
mechanism, if let run to steady state in a closed system 
with radiative losses, will lead to a particle DF that is 
peaked rather than power law. At first blush this seems 
to violate observations which suggest a broader, probably 
power-law DF. This is, however, not a problem: the mag- 
netic field is almost certainly inhomogeneous (c./. Eilek 
1999 for arguments on field filamentation) . Synchrotron 
emission from a peaked peaked particle DF in a nonuni- 
form field will lead to a broad photon spectrum. Eilek 
& Arendt (1996) calculated this for a choice of particle 
and field distributions. Figure 4 shows some examples of 
broad synchrotron spectra produced by a narrow particle 
DF. In addition, the particle DF itself may be field depen- 
dent. The critical energy pc will reflect the net held seen 
by the particle in its lifetime, so that an inhomogeneous 
magnetic field will lead to a range of pc values. This will 
also lead to a broadening of the final photon spectrum. 

Finally, we might ask, how common are radio haloes in 
clusters? We suspect that all clusters contain magnetized 
plasma, and also have relativistic particles injected from 
galaxies in some way or another. The simple conclusion 
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Fig. 4. Examples of broad synchrotron spectra from nar- 
row particle distributions. The fraction of the source vol- 
ume filled with magnetic field at strength B is denoted by 
g{B). The solid line has f{p) oc 5{p—po), and g{B) oc B~^ 
for Bi < B < B2. The dotted hne has f{p) a p^e-^'P", 
and g{B) oc B~^ for B < B2 {Bi, B2 are arbitrarily cho- 
sen cutoffs). The dashed line has f{p) oc S{p — Po), and 
g{B) oc Bi/^e-B/Bo. Taken from Eilek & Arcndt (1996). 
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is that the ICM in all clusters is a source of diffuse syn- 
chrotron emission. However, it is becoming clear that both 
particle acceleration and magnetic field strength (c/. Eilek 
1999) are sensitive to the level, and nature, of turbulence 
in the fCM. The bolomctric synchrotron cmissivity is, of 
course, 



jsy ^ -necaril) B oc u^B 



(12) 



where (7) is the mean energy of the relavistic electrons, 
rie is their number density, and Ue = rie (7) is their energy 
density. Each of these quantities, as well as the magnetic 
field, is determined in part by the underlying turbulence. 
Thus, synchrotron emission is merely an indirect tracer of 
the turbulence, and one with a very sensitive gain factor. 
Perhaps both the turbulence and its gain factor varies 
from cluster to cluster, and we have only detected the 
brightest haloes and relics up to now. 
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